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Abstract 

The important parameter (Iq of the heavy quark expansion is analyzed including 
perturbative and power corrections. It is found that /iQ(2GeV) is known with a 
few percent accuracy. The perturbative corrections are computed and found small. 
A nonperturbative relation is suggested which allows to control the power correc- 
tions. We conclude that yU^(l GeV) = 0.35 to;o2 GeV^. The two-loop expression for 
the effective "M£"' radiation coupling a^/^^^cu) is given which improves reliability 
of the perturbative evolution of fiQ towards the low momentum scale. On the non- 
perturbative side, we advocate the utility of combining the heavy quark expansion 
with expanding around the "BPS" -type approximation for the meson wavefunction, 
which implies relations /i^ ~ fiQ and —p^s ~ Pz) ^^^^ similar ones for the 
nonlocal correlators. 
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The heavy quark expansion allows to quantify the effects of nonperturbative 
physics in beauty decays, often in a model-independent way starting from the first 
principles of QCD. The most informative predictions are obtained for observables 
where the Operator Product Expansion (OPE) applies, like the inclusive decay 
widths. The leading nonperturbative effects are described by the two heavy quark 
expectation values /x^ and fiQ of the kinetic and chromomagnetic operators, re- 
spectively, entering at the Aqq^^/^I ^^^^^ 0- It turns out, in particular, that the 
extractions of Vcb are sensitive to them. 

The value of (Iq can be determined from the mass splitting between B* and B 
mesons. The kinetic expectation value fi"^ is a priori more uncertain. The set of 
heavy quark sum rules, however, significantly restricts /x^ in terms of fiQ [0, ^. In 
particular, fi^ — fiQ is positive, yet unlikely to exceed 0.2 GeV^. 

In this paper we analyze the value of fiQ. Throughout the paper it it refers to 
the infinite mass limit. In QCD the heavy quark operators depend on the renormal- 
ization point. We consistently take this into account using the complete definition 
of these operators in the framework of the quantum field theory suggested earlier 
To this end the perturbative correction to the hyperfine splitting is computed. 

We also address the power corrections. A nonperturbative relation is suggested 
which looks well satisfied in QCD. It allows one to control a number of higher-order 
effects. We conclude that the associated uncertainty in fi^ is rather small. 

We propose and briefly discuss another theoretical tool based on the smallness 
of the difference /i^ — /i^, which can shed light on the structure of the higher-order 
power corrections and improve the accuracy of certain applications of the heavy 
quark expansion to charm. 



1 Perturbative effects 

The value of fiQ^fi) is extracted in practice using the relation = ^mb(MB*—MB)- 
There are perturbative corrections to this relation depending logarithmically on 
n/rrib. At finite rrib there are also nonperturbative corrections suppressed by powers 
of l/mfo. 

The bare heavy quark operator Og = Q-^gsO'f^uG^'^Q = —QgscrBQ is ultraviolet 
divergent in the non-Abelian theory. The normalization point is introduced via 
the upper cutoff in the integral over the antisymmetric small velocity (SV) heavy 
quark structure function W-{e) expressing the sum rule for the matrix elements of 
/ig. In the standard heavy quark notation it takes the following form: 

^= rW4e)e'de = 2 ^ el^M - ^ el\nM\ (1) 

where r's are the P-wave transition amplitudes and Sk are the corresponding exci- 
tation energies (for a review, see 0). The OPE gives the SV structure functions in 
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terms of the zero-recoil matrix elements of the momentum operators -Kj = QiDjQ: 

n 

Since \D^^D]\ = —igsGji = itjikQsBk holds, the sum rule (0) is transparent 0|: 
Qtt^Q = Efc TTfcVTfc , and QqsBiQ = ieijkiTjTik- 

In QED the integral in Eq. (|lD converges and defines the magnetic field strength 
ei?em(0) at the position of the static center. The magnetic spin interaction of an 
elementary heavy fermion is given precisely by this expectation value times the Dirac 
anomalous moment, eB^j^a ^(l + ^ + In a non-Abelian theory like QCD the 

integral diverges in the ultraviolet, and the expectation value of (7s-Bchr(0) depends 
on the normalization point. Let us note that the adopted definition of the operator 
corresponds to the usual scheme with the two covariant derivatives taken at different 
points and connected by the P-exponent. The displacement lies on the (Euclidean) 
time axis and its magnitude is governed by 1//^. More precisely. 

The heavy quark Hamiltonian has the well known form 

nQ = mQ-g,A, + ^ + ca^^^^^ + o(^] . (3) 
2mQ 2mQ \mf^J 

To compute perturbatively the Wilson coefficient ccifJ') we therefore can consider 
the zero- velocity heavy quark transition amplitude Tij{u) mediated by the currents 
QiDiQ and QiDjQ, Fig. 1 (the corresponding OPE formalism is discussed in detail, 
e.g. in Ref. [^). To select the matrix element of the chromomagnetic operator we 
evaluate it on a heavy quark state including scattering of an additional (transverse) 
gluon, and look for the component antisymmetric in This appears in the linear 
in the gluon momentum q approximation. For simplicity, we take qo = 0, and assume 
the initial quark is at rest. We also do not show explicitly the heavy quark spin 
indices (as if Q is a scalar). 




Figure 1: The heavy quark forward scattering amplitude as a function of energy to. The 
solid blocks denote the momentum operator QiDQ. 
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The tree level 0{a^) expression is obvious: 



Tijiuj) 



-Lu + ie 



gsgj^iiQ^Q 



(4) 



where I and c are the gluon polarization and color indices, respectively. The tree 
gluon QCD vertex is gs^pQ7^.Y'pQ = 9s^Pq {^^^^^ ' ^) H^Q- Comparing its 



nonrelativistic expansion ipQjitpQ ~ {Pi+P2)iVqVq + if^ijkVq^^j^QQk with T^^ 
yields c^q"^ = 1. 

To account for the strong interaction corrections we compute Tij{uj) perturba- 
tively assuming —u ^ Aqcd- The same corrections are computed for the QQg 
vertex in QCD projected onto the magnetic spin structure. The vertex yields the 
one-gluon matrix element of the effective heavy quark Lagrangian. The difference 
between the two determines the coefficient cg- It is ultraviolet (UV) finite, as well as 
infrared (IR) finite even at g ^ 0. The latter limit significantly simplifies computa- 
tions allowing for directly expanding the Feynman integrands over q. The resulting 
integrals are saturated in the domain of momenta between u and mg. The diagrams 
one has to compute in the static theory are shown in Figs. 2. 







Q 




Figure 2: Examples of one-loop diagrams for the one-gluon matrix element of Tij{Lv). 
The remaining diagrams vanish for the chosen kinematics, or upon antisymmetrization 
over The gluon wavefunction renormalization is omitted. 

The individual diagrams, however, can be ultraviolet and infrared divergent. 
Since the whole integral for cg is well behaved, one can use any regularization in the 
infrared and the ultraviolet for computing separate diagrams; the only requirement 
is that it must be consistently the same. For example, one can compute cg in 
dimension D = 4: + 2e; then at e > there are no IR singularity and the limit 
g — >■ is straightforward.^] On the other hand the point e = is regular for cq- 
However, since there are separate diagrams which diverge both in the UV and IR, 
dimensional regularization is not advantageous. Instead, we cut the integrals in the 

^At D = A the chromomagnetic moment in QCD has mfrared divergence ^ Inl/g^. 
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UV at /c^ = A^, and introduce the IR mass regulator in the gluon propagator ^tzj2- 
The latter allows us to use the limit —>-0 even at D = 4. Cancellation of the terms 
dependent on and provides a useful cross-check. 

The computation of the one-loop matrix element of Tij in the effective static 
theory results in the Feynman gauge in 



iTij-Tji){uj) = {SnQj - SjiQi) gsQ^Q 




(5) 



The QCD vertex takes the form 



7/. 



-Ca 



167r2 



21n-- 



7/. + 



Cf + 



-Ca 



In 



mi 
A2 



2mQ 



(6) 

where Cp = 4/3 and Ca = Nc- We have omitted from both expressions the dia- 
grams renormalizing the external gluon propagator, since they are the same in both 
cases. For magnetic structure the difference in the Abelian part amounts, as ex- 
pected, to the Schwinger anomalous term Cp^- The no n- Abelian difference is just 

-C^^ln^. 

We note, however, that with the non- Abelian interaction the two theories would 
have different running strong coupling gs{k'^) below thq if the bare coupling g^^^ and 
the UV cutoff A are taken the same. This is seen by evaluating the "charge" gluon 
vertex Fq in the static theory. The difference originates from the part of the 'Abelian' 
vertex correction proportional to Ca which is not canceled by the renormalization 
of the quark wavefunction - it yields pure In • (The 'non- Abelian' vertex is absent 
from the charge interaction.) The total correction to the 7^ structure in the vertex 
in QCD does not depend on the quark mass, Eq. as it should be to respect gauge 
invariance. However, this holds only provided the UV cutoff is infinitely larger than 
all other mass scales. The renormalization of the static quark interaction differs by 
a constant since here the UV cutoff is much lower than mg. This means that the 
static quark gauge interaction requires a different counterterm. Alternatively, it can 
be expressed by saying that for static quarks Agtat niust be taken different from A in 
full QCD. As follows from Eq. (|^), at one loop one has Astat = A/e2 in the Ca term, 
and ^ Ca In ^ gets replaced by ^ C/ 



2lT 



(in ^ - 2 



To calculate the matrix element of ^j.q{h) we need to integrate ^^^^^..^ij 



lmTij{u). The 
even in perturbation theory. However, 



literal expression (||) is not valid at |ti;| ^ \q 
the perturbative T{uj) has the proper analytic properties to any order. We then can 
represent the sum in Eq. (|l|) as an integral over the contours in the complex u plane 
stretched away from small u, see Fig. 3 where the perturbative matrix element is 
given by Eq. (^. The integral is simple: 



{Qg\gsGij{fi)\Q) = g,{6iiqj 



1 - 



-Ca 



167r2 



4 In 



A 



(7) 
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We thus get the final result 



/i 



CG(/i) = 1 + Cf-^ + Catt In — + 2 
2n 2% \ niQ ^ 



(8) 




Figure 3: The complex plane of energy uj. Thick line at cj > shows the cut of T(uj). 
The integral of Imr(a;) can be taken over the circle =/i. 

It is convenient to absorb the constant term 2Ca in the non-Abelian part into 
the argument of the logarithm, it i.e. use (In + 2 ) = In Moreover, the usual 
Abelian part ^Cp depends on the normalization convention used for the heavy 
quark mass. The standard Schwinger coefficient in Eq. refers to the pole mass, 
a choice disfavored in QCD. Using instead the running 'kinetic' mass mQ{jl) ^ 
we find that for the chromomagnetic term in the Hamiltonian (|^) 



caM _ l + gAtlngg 
2mQ 2mQ(^) 

Alternatively, using the MS mass this expression becomes 

c' (lA l + CAT^ln^ 

'^G\f^) _ 27r mq _ 



(9) 



2mQ 2mQ{mQ/e3) 



(10) 



however, the MS mass of a heavy quark loses physical significance at the normaliza- 
tion scales below mg 0. 



2 Power corrections 

The perturbative relation 

CG{fi)fiUf^) = ^mi,-{MB*-MB) (11) 

holds only for asymptotically heavy b where power-suppressed effects die out. The 
1 / ml corrections to the hadron masses are given by two local heavy quark operators 
and four nonlocal correlators also describing the l/rriQ corrections to the hadrons' 
wavefunctions. For Mb*— Mb one has 
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The expectation values of the convection current (or spin-orbital) operator p^^ = 
^^{B\cr-Exn\B) would vanish to the extent that B could be described as purely 
a two-body system like in nonrelativistic approximation. We know, however, that 
in actual QCD the heavy quark bound state is rather relativistic. This is quantified 
by the diff'erence between the transition amplitudes T3/2 and T1/2 and between the 
masses of the | and | P-wave states, which are the same in nonrelativistic systems 
0. In particular, for the first three moments we have 



0.7 (13) 
0.7 (14) 
0.8 (15) 
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with A ~ 700 MeV and S ~ 250 MeV. The normalization scale dependent A, 
and are taken at the scale around 1 GeV. 

The magnitude of p\g can then be estimated using the next spin sum rule 0: 

- 0ls ~ /i^Phadr ^ 0.15 to 0.2GeV=^ , (16) 

where //hadr ~ 500 MeV is a characteristic mass scale for the P-wave excitations. 
The nonlocal correlators are saturated by the transitions into the = |^ for fi^Q 
and into = 7^ ■, J'^ = |^ for (T'j^ "radial" excitations; they are largely unknown. 

We can estimate the necessary combination of the above spin-triplet D = 3 pa- 
rameters employing the empirical observation that the mass-square splitting between 
vector and pseudoscalar mesons is nearly a constant: 

Ml -Ml^ - ^ Ml, -Mlc:i Ml. - Ml (17) 

which extends even to strange charmed mesons. (A 12% decrease for B fits well 
the expected perturbative renormalization.) It is related to the universal slope 
of the corresponding Regge trajectories, a yet poorly understood nonperturbative 
phenomenon of strong dynamics, perhaps related to a certain simplification in the 
large limit. This universality must be definitely violated for very heavy quarks 
due to hard gluons with momentum scaling with mq. Rather, it can be viewed 
as an inherent property of soft nonperturbative interactions responsible for physics 
around 1 GeV scale. 

The universality implies the relation 



-pi5 + p'G + Pi)^2A/i^^0.5GeV 
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where we have used A ~ MB — mb — 700 Me V. The above nonperturbative param- 
eters are normahzed at the scale around 1 GeV. With the estimate (|TBp we are led 
to an evaluation 

- (pIg + Pa) ~ 0.6 to 0.7 GeV^ . (19) 

The scale of these correlators lies in the expected range if one bears in mind the 
magnitude of other parameters A, yU^, fi^ all being given by a mass 600 to 700 MeV 
to the corresponding power, though possibly on the upper side. 

The sign of p^Q and p\ is not known a priori. We note, however, that the so far 
observed nonperturbative phenomena in the heavy mesons fit well the approximation 
that the ground state B has nearly the "lowest Landau level" wavefunction, or is 
the BPS-saturated state. For instance, /i^— /i^ is noticeably lower than Pq. If 
this saturation was actually the case and fi'^ — fiQ = 0, one would necessary have 
(TTT \B) =0, meaning that the asymptotic wavefunction of the light cloud annihilates 
two certain linear combinations of the total momentum operators, e.g. Vz and Vx—iVy 
for the state jz = \- Signifying vanishing of all ri/2, this would also entail a series of 
relations which include vanishing of all the B meson correlators involving powers of 
avr . Say, 

-pIs = pI^ pIg = -'^pL, pa + pIg = ~pL~pI (20) 

would hold. The first relation agrees with the estimate (|r^) (note that — > Pd)- 
The correlators p^^ and p|, on the other hand, are positive, which shows consistency 
of the estimates. 

Adopting relation (|18|) as the guideline, we end up with 



pl = lca\p;mb)-{M],.-Ml) (21) 

through order l/rrif,- The effect of the 1/mf, corrections in Eq. (jl^) amounts to 
the apparent decrease of p^ by about 15%, recovered in the analysis above. The 
magnitude of the shift is moderate and fits well the expectations 

3 Numerical estimates 

Applying equations (|]), (|12]) and (|T8D we arrive at the evaluation 

p^(2.5 GeV) ~ 0.312 GeV^ , (22) 

where the value of mft(l GeV) = (4.57 ±0.05) GeV was used and its evolution to the 
scale 1.5 GeV ~ m^/S can be safely done using 

dm,(p) ^16 ^ 4 p \ a(f\p) ^^3) 



dp V 3 3 nib J TT 
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the effective dipole radiation strong coupling ai'^^ is known to two loops 0. The 
uncertainty in the mass affects iiq here only at a percent level. Likewise, the per- 
turbative effects are included up to the second loop; their uncertainty can hardly be 
significant. 

The possibly largest uncertainty comes from the precise value of the D = 6 
operators. It has been mentioned above that the 12% decrease in M^*—M'^ compared 
to M^f—M"^ is well described by the perturbative renormalization amounting at one 
loop to a factor of about 0.87. However, the two-loop result by Czarnecki and Grozin 

further suppresses Cq /cq down to about 0.8. Taken literally, this would imply 
even some overshooting, so that the power corrections governed by —{p^g~^ Pa~ Pls) 
have to make up for the extra suppression. This would imply the combination 
to further exceed the estimate ([T8|) by a factor as large as 1.5 to 2, being next to 
fatal for the expansions in l/rric. However, we note that the major perturbative 
factor in the extra NLO enhancement of Mr,*— Mo comes from the increase in the 
renormalization of the chromomagnetic moment of the charm quark due to growing 
Us at the charm scale. Introducing the proper Wilsonian cutoff around 1 GeV would 
effectively stop the increase of this short-distance coefficient for charm, and eliminate 
the apparent contradiction with the observed behavior without invoking a too large 
scale for higher-order power corrections. 

To remain on the safe side we allow for an additional factor 0.75 to 1.4 in Eq. (0). 
Then we arrive at 

/i^(2.5 GeV) = (0.30 to 0.33) GeV^ . (24) 

The scale 2.5 GeV would be high enough to trust the perturbative expansion. In 
practical apphcations we need, however, to evaluate /i^ at the lower scale around 

1 GeV. The perturbative evolution of Pq obeys 

^d^|^o_^^fi:^^.,^,_ (25) 

where a^™*^), in analogy with the coupling in Eq. (p3D can be called "ME"' -radiation 
one, or "ME-coupling" . To first order it is the usual QCD a^. The complete two- 
loop calculation presented in paper combined with the one-loop result Eq. (H) 
is sufficient to determine ^^'"^•'(ti;) to order a^, using renormalization invariance of 
observables. This yields 

Ca'^^ = 0V„ - 2C4JiC. - fn,) (^) ' + O (af ) , (26) 

where 7^ is the anomalous dimension of Ref. 0, Eq. (15) and nf is the number of 
light flavors. Therefore, we have 



1 25 ^ 
1 - 2i'-'A — 



a, 

TT 



Ca = iV, = 3 . (27) 
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As expected, the two-loop anomalous dimension of the physically defined chromo- 
magnetic operator differs from the MS one already at order (its conformal part 
is still the same to two loops). 

The evolution equation Eq. (E^l) is immediately solved: 



Ca 

TT I Jh_ 

+ 2/3„ 



f28) 



where (3q and (3i are usual QCD one- and two-loop coefficients in the /3-function, and 
the evolution of the ME coupling to two loops is given by the same renormalization 
group equation. 

The evolution equation (|25|) viewed perturbatively suggests that yUG(yu) is en- 
hanced toward lower normalization scale. Taking it at face value we get an enhance- 
ment factor 1.11. This is larger than the literal one- loop difference coming from 
Eq. (^) due to the growing strong coupling which is accounted for in the evolution 
equation. One clearly cannot go with the perturbative evolution too low in /x. For 
instance, aX ji^Q the value of 11% goes to zero, constituting apparently only about 
0.2 GeV^ at /i = 0.5GeV 0. Nevertheless, it is justified to trust the above mod- 
erate enhancement. The physically defined coupling is smaller than and 
remains suppressed protecting perturbative corrections from blowing up. This is in 
contrast with the standard MS anomalous dimension which receives large positive 
two-loop correction. 

Combining the above factors, we arrive at 

/i^(lGeV) = | (m|.-M|)-[0.94±0.025 ^1^35 ]• (29) 

ZL \ / pert power 

Therefore, we conclude that the value of GeV) most probably amounts to 

/i^(lGeV) = 0.35^°:°^GeV2; (30) 

the lower values can appear in an unfavorable scenario where power corrections in 
charm go out of theoretical control. 



4 Discussion and outlook 

The presented analysis allows us to determine the consistently defined chromomag- 
netic value /i^ normalized at the scale around 2 GeV with minimal theoretical un- 
certainty. We conclude that the value /i^(l GeV) ^0.4GeV^ we have used so far 
was reasonably accurate, yet probably about 10% larger than it is in reality. In 
principle, as low a value as 0.30 GeV^ cannot be rigorously excluded at present, for 
the price of endangering the l/rric expansions. The uncertainty in the perturbative 
effects can be further decreased carrying out the same program to two loops. A 
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complementary information on the D = 3 nonperturbative parameters in the hadron 
mass expansion would be helpful to get more confidence in the evaluation of the 
power corrections and to shrink the error bars down to a few percent level. 

We found the essentially non-Abelian effective "M£"' -coupling which does not 
have a counterpart in electrodynamics, either classical or quantum. It determines, 
for example, the fraction of the total momentum of the QCD degrees of freedom 
carried by the hard (short- distance) components of the light cloud; it is power sup- 
pressed, cx [0. Physics behind this phenomenon will be discussed elsewhere. 

At the loop level the gluon self-interaction suppresses the effective "M£"' -coupling 
compared to the traditional MS coupling. The conformal part of it turns out similar 
to the one in the dipole radiation coupling 0, although the non-Abelian suppres- 
sion is approximately twice larger then in the latter. If iVc — vs. ( if — t| ) A^c — • The 



24 c ^ '^^ y Q 12 

similarity further supports the a priori expected advantage of using the radiation 
coupling as an effective perturbative expansion parameter. Following the OPE-based 
line of reasoning of Ref. we find that the l/u nonperturbative component in this 
coupling at large u is associated with the LS operator and is proportional to its 
anomalous dimension 7^5. Implementing the approach of Ref. one can show 
that this anomalous dimension exactly coincides with that of the chromomagnetic 
operator. Q Then we conclude that in B mesons, for instance, at large u 

S^^aT^Hu) ^ -pUu) ^0.5GeV ^^^^ 



a 



This estimate is consistent with the possibility to perturbatively evolve /i^ down to 
a 1 GeV scale. 

The precise expectation value of the kinetic operator /x^ is quite critical in a 
number of applications. Heavy quark sum rules ensure that the inequality > 
holds at arbitrary normalization point 0, and /i^ almost certainly must 
lie in the interval 0.4 GeV^ < /i2(lGeV^) < 0.55 GeVl Moreover, using the spin 
sum rules [0 one has — = 3e^- — 0.75) with 0.5GeV^e</i, and 

we expect £~0.5 GeV for the usual choice of /i = 1 GeV. (One should keep in mind 
that the value of the kinetic energy of actual h quark in B meson can be reduced 



due to often discarded, but probably noticeable l/rrib effects p, pA|.) 

We anticipate that combining the heavy quark expansion with an approximation 
assuming fi^ — fiQ <^ fiQ manifesting the proximity to the "BPS" regime for the 
ground state, can be useful in guiding us through understanding pattern, or even 
physics of higher-order power corrections in B and D mesons. Expanding around 
this approximation provides a new and effective nonperturbative parameter small 
enough to isolate a number of potentially large power corrections. For example, the 



■^The anomalous dimensions of cq p^^ and Cq p\ vanish by the same token, and their mixing 
into pls is easily fixed. I am grateful to M. Eides and A. Vainshtein for discussing this point. 
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usual spin-averaged heavy meson mass expansion replaced by 



(32) 

is possibly more stable in respect to higher orders being governed by the suppressed 




higher-dimension expectation values scaling like powers of y/x^ — /x^. As has been 
pointed out , the uncertainty in the precise value of mb—mc is currently the main 



limiting factor in extracting \Vcb\- The validity of the approximation can be exper- 
imentally cross checked by carefully analyzing the semileptonic b ^ c transitions 
into excited states, in particular to the ^ P-waves. 

In the BPS regime, a number of other nonperturbative effects become more 
tractable. Say, the B^D formfactor at zero recoil does not have I/'ht-q corrections 
at all. The corrections still are present and significant in the B —>■ D* zero recoil 
formfactor -F(O) , but their structure simplifies: 

f(0)=d-(l+X)||-o(4) (33) 
with X Pg given by a sum of two other positive nonlocal correlators p^^ + p| 



4Mb 



{B\iT{b7r^b{x),b7r^b{0)\By+ -{B\iT{bajBkb{x),bajBkb{0)\By 

3 



similar to pl^ and p| (see Eqs. (28) of Ref. 0), but containing the extra factor 
of il^ol in the integrand, or l/(i?„ — i?o) in the language of perturbation theory 
in quantum mechanics. We expect the approximation p^^ + p| ^ (p^^ + p|)/5rd 
(actually, a rigorous upper bound 0) to be reasonably accurate, with ~600 MeV 
the energy of the first j^= radial excitation of the ground state. Taken literally 
this would suggest x to be quite large, around 2. 

The limit p^ — p^ = means that the heavy quark wavefunction minimizes the 
momentum square operator in a given chromomagnetic field. This happens for the 
lowest Landau level which is an example of a "BPS-saturated" state. It is worth 
noting that the newer relativistic quark models of heavy hadrons [0 properly imple- 



menting Lorentz transformations yield a good approximation to this limit. Comple- 
mentary to this, the Block & Shifman QCD sum rule analysis of the IW function 



strongly supports this by virtue of the spin sum rules. In usual quantum mechan- 
ical systems of electrons the BPS saturation is realized applying strong magnetic 
field. In mesons the chromomagnetic field is a priori of the same order Aqqj-, as the 
chromoelectric field, and is far from classical. In B mesons the approximate BPS 

^It is interesting to recall that the limit of a strong magnetic field in a quantum mechanical 
system of charged particles yields a physical realization of noncommuting space coordinates [|l5j. 
Actual B mesons may well share some of their properties. 
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saturation would rather be dictated by a very special internal structure of the light 
cloud, leading to the strong correlation between the spin and momentum operator. 
Such a correlation vanishes in a nonrelativistic system, and can be realized only 
in a deeply relativistic regime. It does not look probable that such a property, if 
confirmed experimentally, is purely accidental. Perhaps, it is related to a certain 
large parameter, like the number of space dimensions or the number of colors. It is 
intriguing to study such possible connections. 

Our analysis of fiQ incorporating powerful constraints from a series of the heavy 
quark sum rules gives an indication that the higher-dimension expectation values are 
quite significant and probably lie at the higher end of the existing estimates. Like- 
wise, the pattern of the excited heavy quark states seems to be at some variance 
with the routinely employed assumptions inherited from obsolete models. There are 
fresh ideas of how to improve the knowledge in less vulnerable ways. These will be 
addressed in forthcoming publications. 
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